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Abstract — Excitation of radial pulsations in red supergiants of Magellanic Clouds is in- 
vestigated using the stellar evolution calculations and the self-consistent solution of the equa- 
tions of radiation hydrodynamics and turbulent convection. The stars with initial masses 
6M Q < M Z ams < 28M and the initial chemical composition X = 0.7, 0.004 < Z < 0.008 



are shown to be unstable against fundamental mode oscillations with periods from 17 to 1200 



days as they become helium burning red supergiants. The period-luminosity relation slightly 
depends on the mass loss rate varying with a factor of three, whereas its dependence on the 
metal abundance is given by 5M bol = 0.895 log Z. In comparison with galactic red supergiants 
the low metal abundances in red supergiants of Magellanic Clouds are responsible for their 
higher effective temperatures and substantially narrower ranges of evolutionary radius change 
during helium burning . Therefore on the period-mass diagram the red supergiants of Magel- 
lanic Clouds are located within the strip with width of 5\ogM 0.09, so that the uncertainty 
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of mass evaluation of the red supergiant with the known pulsation period is nearly 25%. 
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Red supergiants belong to brightest stars [L < 3 x 1O 5 L ) and are also known as long- 
period variables (LPV) with cyclic light variations on the time scale from a few dozen to ~ 10 3 
days. The pulsational nature of light variability in red supergiants has been firstly shown by 
Stothers (1969) and later was supported in several following works (Stothers 1972; Li and Gong 
1994; Heger et al. 1997). In the General Catalogue of Variable Stars (Samus et al. 2011) these 
stars are classified as pulsating variables of type SR with semiregular light variations. 

The pulsation hypothesis for the light variability of red supergiants is strongly supported 
by the period-luminosity relation derived by Glass (1979) from photometric observations of 
38 brightest M-type supergiants in the Large and Small Magellanic Clouds (LMC and SMC). 
Since then the number of semiregular late-type variables observed in these stellar systems 
increased by several times (Yang and Jiang 2011, 2012). The period-luminosity relation for 
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red supergiants is of great importance due to the fact that in comparison with Cepheids they 
are more appropriate for extragalactic distance calibration (Pierce et al. 2000; Jurcevic et al. 
2000). 

However in contrast to Cepheids with well known mechanism of pulsation instability the 
red supergiants still remain scarcely studied pulsating variable stars. Theoretical studies of 
pulsational instability in red supergiants were presented only in two papers by Li and Gong 
(1994) and Guo and Li (2002) devoted to the linear nonadiabatic analysis of radial oscilla- 
tions. The authors concluded that the observed variability of red supergiants is due to radial 
oscillations in either the fundamental mode or the first overtone and pulsations are driven by 
the K-mechanism in the helium ionizing zones. Unfortunately, the linear analysis employed in 
these studies is based on the simplified treatment of convection and does not take into account 
effects of turbulence (Li 1992a, 1992b). In red supergiants the turbulence plays an important 
role because of the large extension of the outer convective zone (as large as 70% of the stellar 
mass) as well as due to substantial contribution of the turbulent pressure which is as large 
as one third of the total thermodynamic pressure (i.e. the sum of the gaseous and radiation 
pressure). 

Hydrodynamic calculations of nonlinear radial oscillations in red supergiants based on the 
self-consistent solution of the equations of radiation hydrodynamics and turbulent convection 
were done by Fadeyev (2012). This study deals with stars with initial masses 8M < M Z ams < 
20M Q and relative mass abundances of hydrogen and elements heavier than helium (metals) 
X = 0.7 and Z = 0.02, respectively. However, this composition is typical for Population I 
stars of our Galaxy, whereas in stars of LMC and SMC the metal abundancies are significantly 
lower: 0.003 < Z < 0.008 (Peimbert and Torres-Peimbert 1974, 1976). In the present paper 
we show that lower metal abundances in red supergiants of LMC and SMC lead to effective 
temperatures higher by 300 K < AT c g- < 500 K in comparison with galactic red supergiants. In 
late-type stars (T eff ps 3500 K) such a difference in effective temperatures leads to significant 
changes in both stability of the star against radial oscillations and the pulsation period. 

The goal of the present work is to study the pulsational properties of red supergiants in LMC 
and SMC with methods of hydrodynamic computations of nonlinear radial stellar oscillations. 
We consider the stars with initial masses 6M < Mzams < 28M Q and relative mass abundances 
of hydrogen and metals on the zero age main sequence (ZAMS) X = 0.7 and Z = 0.004, 0.008, 
respectively. 
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THE METHOD OF COMPUTATIONS 
STELLAR EVOLUTION 

The study of self-exciting stellar oscillations is the Cauchy problem for equations of radia- 
tion hydrodynamics and turbulent convection that describe the spherically-symmetric motion 
of the self-gravitating gas. Initial conditions in such a problem are the spatial distributions of 
physical variables corresponding to the hydrostatic equilibrium. In the present work the initial 
conditions are obtained from the stellar evolution calculations. Description of the methods 
employed are given in our previous papers (Fadeyev 2007, 2010). 

In calculations of stellar evolution we used the steady-state local convection model (Bohm- 
Vitense 1958) with the ratio of mixing length to pressure scale height t/H-p = 1.6. It was also 
assumed that the size of the convective core increases due to convective overshooting by O.lifp. 

Stars with initial masses 6M < M Z ams < 15M occupy the red supergiant domain only 
during the initial phase of core helium burning and later they leave this domain as helium is 
exhausted in the convective core. For example, the star M Z ams = 1OM , Z = 0.008 leaves the 
red supergiant domain when the central helium abundance is X( 4 He) c < 0.4. In such a case 
the evolutionary track loops the Hertzsprung-Russel (HR) diagram and crosses the Cepheid 
instability strip. More massive stars (15M Q < M ZA ms < 28M ) remain red supergiants during 
the whole stage of thermonuclear core helium burning and the shape of their evolutionary tracks 
in the HR diagram depends on the mass loss rate M. 

Various empirical formulae for M were recently discussed by Mauron and Josselin (2011) 
who concluded that the best agreement with observations of red supergiants in Magellanic 
Clouds is obtained with the formula by de Jager et al. (1988). In the present study the mass 
loss rate of evolving red supergiants was assumed to be M = /M NJ , where M NJ is the mass 
loss rate evaluated by formula of Nieuwenhuijzen and de Jager (1998) which only insignificantly 
differs from that by de Jager al. (1988). The factor / is introduced to estimate the dependence 
of results of hydrodynamic computations on uncertainties in mass loss rates and in stellar 
evolution calculations was fixed in the range 0.3 < / < 1. 

The role of the mass loss rate in the evolutionary track of the red supergiant with initial 
mass Mzams = 2OM and metal abundances Z = 0.004, 0.008 and 0.02 is illustrated in Fig. I. 
Of most interest are parts of the track that correspond to the helium core burning and in Fig. 
1 they are shown by solid lines. Dottedd lines correspond to the evolutionary stage when the 
energy source in the stellar center is gravitational contraction of the core and movement in the 
HR diagram is several orders of magnitude faster than during thermonuclear helium burning. 
We do not discuss the stage of gravitational core contraction because of negligible probability 
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to observe such stars. 



NONLINEAR STELLAR PULSATIONS 

In solution of the equations of hydrodynamics we assumed that the flux of enthalpy F c and 
the flux of the turbulent energy of convective elements F t are defined according to Kuhfufi(1986). 
In contrast to previous works of the author (Fadeyev 2011, 2012) the equation of motion and 
the energy equation are written in the following form: 

d 2 r GM r t 2 d ,„ „. • 
dE dV 1 

aF + p ar = -p V -< F ' + F <>- e " < 2 > 

^ + P t ^ = -Iv.F t+ £ l + e ,. (3) 
at at p 

The equation of motion ([1]) describes the change of the radius r and velocity U = dr/dt 
of the mass zone with Lagrangean coordinate M r due to gravity (here G is the gravitational 
constant) and gradients of the total thermodynamic pressure P (the gas pressure and radiation 
pressure) and the turbulent pressure P t . The last term in the right-hand-side of equation (JT]) 
accounts for the momentum transfer between the gas flow and turbulent elements. 

Equations (j2J) and ([3]) describe the change of the specific internal energy of gas E (the sum 
of the translational energy of gas particles, excitation and ionization energy of atoms, energy 
of radiation) and the specific turbulent kinetic energy E t . Here p = l/V is the gas density, V 
is the specific volume, F r is the radiation flux calculated in approximation of radiation heat 
conduction. The coupling term C t is defined as 

C t = S t - D t - 2) rad , (4) 

where S t is the rate of specific turbulent energy generation due to buoyancy forces, T> t and 
D ra d are the rates of turbulent energy dissipation due to molecular viscosity and radiation, 
respectively. Formulae for these quantities are given by Wuchterl and Feuchtinger (1998). 

The turbulence is assumed to be isotropic, so that viscous stresses are determined by the 
Reynolds stress tensor (see, e.g., Pope (2000)). In spherical geometry the rates of momentum 
transfer lt t and energy transfer £ t between gas flows and turbulent elements are scalar trace-free 
parts of the Reynolds tensor (Wuchterl and Feuchtinger 1998). 
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RESULTS OF HYDRODYNAMIC CALCULATIONS 



PULSATION AL INSTABILITY DOMAIN 



Some models of evolutionary sequences corresponding to the thermonuclear core helium 
burning were used as initial conditions in the Cauchy problem for equations (JTj) (J3J) . The role 
of small perturbations was played by errors arising in interpolation of the evolutionary model 
computed by the Henyey method to the hydrodynamical model with Lagrangean mass intervals 
increasing inward as geometric progression. To diminish interpolation errors the evolutionary 
calculations were done with the number of mass zones ranging from 5 x 10 3 to 10 4 , whereas in 
hydrodynamical models the number of Lagrangean intervals was N fa 10 3 . 

Integration of the equations of hydrodynamics ([I]) - ([3]) with respect to time t was accom- 
panied by evaluation of the kinetic energy 



where Uj is the gas flow velocity in the j—th Lagrangean zone, AMj_i/ 2 is the mass interval 
between (j — l)-th and j-th zones. The value j = 1 corresponds to the inner boundary of the 
hydrodynamical model where U\ — and the luminosity satifies the condition dLi/dt = 0. 

Exponential decrease of the averaged over cycle kinetic energy E-^{t) shows that the star is 
stable against radial pulsations. For the criterion of pulsational stability we used the condition 
that the averaged over cycle kinetic energy decreases by at least two orders of magnitude in 
comparison with kinetic energy of the initial perturbation. 

In the case of the growth of kinetic energy the hydrodynamic calculations were carried 
out untill the limit cycle is attained and the averaged over cycle kinetic energy becomes 
independent of time. This condition is fulfilled for models with initial masses Mzams — 15M 
due to small nonlinear effects. In more massive stars the amplitude of limit cycle oscillations 
varies from cycle to cycle and the condition of the constant cycle-averaged kinetic energy is 
fulfilled only approximately within sufficiently large time intervals. 

The mean period of radial pulsations LT was evaluated from the discrete Fourier transform 
of the temporal dependence E-^{t) for both pulsationally stable and pulsationally unstable 
models. Within the initial time interval with linear change of lni^Kmax the instability growth 
rate r\ = lid In E^ max /dt was evaluated, where -EKmax is the maximum kinetic energy attained 
within one oscillation cycle. For all hydrodynamical models with positive growth rate (rj > 0) 
stellar pulsations were found to be in the fundamental mode. 

The change of the instability against radial oscillations in the evolving red supergiant is 
illustrated by plots of 77 in Fig. 2 for several values of M Z ams- The plots are shown versus the 
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relative mass abundance of helium in the stellar center X( 4 He) c because duration of thermonu- 
clear helium burning depends on the stellar mass. In all models the onset of the core helium 
burning takes place for 0.97 < X( 4 He) c < 0.99. 

The growth rate of pulsational instability depends on many factors but the most important 
is the size of the outer convective zone. Luminosity decrease of red supergiants 15M Q < 
Mzams < 20M o during the initial helium burning (see Fig. 1) is accompanied by shallowing 
of the outer convective zone. For X( 4 He) c ~ 0.4 the mass of the convective zone reaches its 
minimum (nearly a half of the stellar mass) and as is seen from plots in Fig. 2 for M ZA ms = 18M 
and Mzams = 2OM the pulsational instability growth rate becomes highest. Further helium 
exhaustion is accompanied by increase of the both luminosity of the red supergiant and mass 
of the outer convection zone, whereas the instability growth rate r\ decreases. 

Typical plots of the gas flow velocity at the upper boundary U and the bolometric light 
Mboi of the red supergiant with nearly maximum growth rate r) are shown in Fig. 3. The 
initial stellar mass, the factor of the mass loss rate and the central helium abundance are 
Mzams — 2OM , / = 0.3 and X( 4 He) c = 0.39, respectively. For the sake of comparison the 
plots of U and M^ Q \ for the model Mzams — 9M , X( 4 He) c = 0.36 are shown in Fig. 4. 

Instability against radial oscillations enhances with increasing Mzams due to higher stellar 
luminosity and increasing nonadiabaticity of stellar pulsations. This is illustrated in Fig. 2 by 
the plot for the model sequence M Z ams — 25M . Here however one should be noted that for 
7] > 0.2 estimates of this quantity become incorrect because of the small number of maxima of 
the kinetic energy during the stage of instability growth. 

PERIOD-LUMINOSITY RELATION 

The change of the pulsation period II of the evolving red supergiant is mainly due to 
evolutionary changes of its radius R because the period of the fundamental mode is proportional 
to the sound travel time from the stellar center to the surface. In the beginning of thermonuclear 
helium burning the luminosity of the red supergiant decreases and the effective temperature 
increases (see Fig. 1) therefore the fundamental mode period diminishes. Following increase of 
luminosity and decrease of effective temperature are accompanied by increase of the period II. 
Such a general property of red supergiants is illustrated in Fig. 5 for three models with initial 
mass Mzams = 2OM and three different metal abundances Z. Each evolutionary sequence of 
hydrodynamical models is represented by filled circles connected by solid lines for the mass loss 
rate factor / = 0.3 and by dotted lines for / = 1. 

The period-luminosity diagram in Fig. 5 reveals remarkable differences between red super- 
giants of our Galaxy and those of Magellanic Clouds. First, for the fixed luminosity galac- 
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tic red supergiants pulsate with longer periods. As noted above this is due to the fact that 
red supergiants with lower metal abundances have higher effective temperatures and smaller 
radii. Second, ranges of evolutionary changes of radial oscillation periods of galactic red su- 
pergiants are appreciably wider. For example, for stars with initial mass Mzams — 2OM and 
Z = 0.02 the pulsation period changes in the range from 510 to 1020 days. With decreas- 
ing Z the range of evolutionary changes of the radial oscillation period shrinks to 20% and is 
380 day < n < 460 day for Z = 0.008 and 370 day < n < 430 day for Z = 0.004. The cause 
of such remarkable differences between red supergiants of our Galaxy and those of Magellanic 
Clouds is due to the ranges of evolutionary changes of the stellar radius. For example, in 
galactic red supergiants with M Z ams = 2OM the evolutionary change of the radius R is about 
one third, whereas in red supergiants of Magellanic Clouds relative evolutionary changes of the 
radius do not exceed 10%. 

The general period-luminosity diagram involving all hydrodynamical models computed in 
the present work is shown in Fig. 6. In this diagram we represent 83 hydrodynamical models for 
Z = 0.008 (6M < M ZAMS < 28M ) and 19 models for Z = 0.004 (1OM < M ZAMS < 24M ). 
In Fig. 6 are plotted also 7 models for Z = 0.02 with initial masses Mzams — 15 and 2OM . 

The role of the mass loss rate M in the mass-luminosity relation becomes perceptible only 
for stars with initial masses Mzams > 2OM . It is seen in Fig. 5 where two model sequences 
for Z = 0.008 are shown for the mass loss rate factor / = 0.3 and / = 1. The mass loss rate 
increases with increasing luminosity (i.e. the initial mass) of the red supergiant. At the same 
time the role of nonlinear effects also increases with luminosity, so that the scatter of values 
of the radial pulsation period II becomes larger. From hydrodynamic computations we found 
that the scatter due to nonlinear effects is roughly comparable with effects of mass loss, that 
is why the models in Fig. 6 are shown for each value of Z independently of the factor /. 

Linear fits for the results of our hydrodynamic calculations in the period-luminosity diagram 
are given by 



where the period II is expressed in days. For Magellanic Clouds (0.004 < Z < 0.008) the 
dependence of the absolute bolometric light Mboi on Z is approximately given by relation 



M bol = -0.337 - 2.820 log n 



M bol = -0.575 - 2.832 log n 



(Z = 0.004) 
(Z = 0.008) 



(6) 



<5M bol = 0.895 log Z, 



(7) 



which slightly differs from <5Mb i = 0.835 log Z obtained by Guo and Li (2002). 
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PERIOD-MASS DIAGRAM 



The most important application of the stellar pulsation theory is the determination of 
the stellar mass with the period-mean density relation. To this end together with the period 
of light variations II one should have the observational estimate of the mean radius of the 
pulsating star R. Unfortunately, for red supergiants such an approach is impossible because of 
too large uncertainties in observational estimates of their effective temperatures. In particular, 
the problem is complicated due to both distorted continuum because of numerous molecular 
absorption bands and circumstellar reddening by dust grains condensing in the stellar wind. 

However there is another approach to establish the relationship between the stellar mass M 
and the pulsation period II without observational estimates of the stellar bolometric magnitude 
and effective temperature. Such a relationship can be determined because the evolutionary 
changes of the luminosity during thermonuclear helium burning proceed within rather narrow 
ranges. For example, for metal abundances 0.004 < Z < 0.008 the evolutionary luminosity 
variations range from 5\ogL = 0.08 for 25M to <51ogL = 0.14 for M Z ams = 6. 

Evolution of the pulsationally unstable red supergiant in the period-mass diagram is shown 
in Fig. 7 for two model sequences with M ZA ms = 2OM and metal abundances Z = 0.004 and 
Z = 0.008, respectively. In the general diagram shown in Fig. 8 all hydrodynamical models 
with metal abundances 0.004 < Z < 0.008 are concentrated within the strip with boundaries 
given by expressions 



Results of stellar evolution and nonlinear radial pulsation calculations show that oscillations 
of red supergiants in the Magellanic Clouds are due to instability of the fundamental mode. 
In the HR diagram the pulsational instability domain encompasses the wide luminosity range 
(2.5 x 10 3 L Q < L < 2.7 x 10 5 ) because both massive stars (M > 10M Q ) and intermediate- 
mass stars (6M < M < 1OM ) are pulsationally unstable. The theoretical period-luminosity 
relation comprises more than one and a half orders of magnitude (17 day < II < 1200 day). 
This result agrees well with classification of the General Catalogue of Variable Stars (Samus et 
al. 2011) where the semiregular variables SR are mentioned as late type giants and supergiants 
with light variation periods from 20 to 2000 days. Observers studying the period-luminosity 
relation of red supergiants in LMC and SMC dealt so far with stars with periods of light 
variation greater than 200 days, whereas the results presented above allow us to conclude that 





where the period II is expressed in days. 



CONCLUSION 
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the period-luminosity relation extends to shorter periods. 

Determination of masses of red supergiants with the period-mean density relation is im- 
possible because of large uncertainties in observational estimates of effective temperatures. 
However, due to enough small evolutionary changes of the stellar radius during the thermonu- 
clear helium burning the masses of red supergiants in LMC and SMC can be evaluated using 
the observational estimate of the period of light variation from the period-mass relation (JSJ). 
The uncertainty of this estimate <51ogM ~ 0.09 (i.e. no more than 25%) is substantially less 
that for galactic red supergiants where the uncertainty of mass determination is around 50% 
(Fadeyev 2012). 

The recent empirical period-luminosity relations for red supergiants in LMC and SMC are 
determined in magnitudes of near-IR photometric bands and uncertain effective temperatures 
of these stars do not allow the bolometric magnitudes to be correctly determined. Furthermore, 
there is considerable discrepancy between photometrical data obtained by different observers 
(see, for example, Fig. 13 in the paper by Yang and Jiang (2011)). That is why the theoretical 
period-luminosity relation in Fig. 6 was not compared with observational data. 

The study was supported by the Basic Research Program of the Russian Academy of 
Sciences "Nonstationary phenomena in the Universe" . 

REFERENCES 

1. E. Bohm-Vitense, Zeitschrift fiir Astrophys. 46, 108 (1958). 

2. Yu.A. Fadeyev, Pis'ma Astron. Zh. 33, 775 (2007) [Astron. Lett 

3. Yu.A. Fadeyev, Pis'ma Astron. Zh. 36, 380 (2010) [Astron. Lett 

4. Yu.A. Fadeyev, Pis'ma Astron. Zh. 37, 440 (2011) [Astron. Lett 

5. Yu.A. Fadeyev, Pis'ma Astron. Zh. 38, 295 (2012) [Astron. Lett 

6. I.S. Glass, MNRAS 186, 317 (1979). 

7. J.H. Guo and Y. Li, Astrophys.J. 565, 559 (2002). 

8. A. Heger, L. Jeannin, N. Langer, et al., Astron. Astrophys. 327, 224 (1997). 

9. C. de Jager, H. Nieuwenhuijzen, and K. A. van der Hucht, Astron. Astrophys. Suppl. 
Ser. 72, 259 (1988). 

10. J.S. Jurcevic, M.J. Pierce, and G.H. Jacoby, MNRAS 313, 868 (2000). 

11. R. Kuhfufi, Astron. Astrophys. 160, 116 (1986). 

9 



. 33, 692 (2007)]. 

. 36, 362 (2010)]. 

. 37, 403 (2011)]. 

. 38, 260 (2012)]. 



12. Y. Li, Astron. Astrophys. 257, 133 (1992a). 

13. Y. Li, Astron. Astrophys. 257, 145 (1992b). 

14. Y. Li and Z.G. Gong, Astron. Astrophys. 289, 449 (1994). 

15. N. Mauron and E. Josselin, Astron. Astrophys. 526, A156 (2011). 

16. H. Nieuwenhuijzen and C. de Jager, Astron. Astrophys. 231, 134 (1990). 

17. M. Peimbert and S. Torres-Peimbert, Astrophys. J. 193, 327 (1974). 

18. M. Peimbert and S. Torres-Peimbert), Astrophys. J. 203, 581 (1976). 

19. M.J. Pierce, J.S. Jurcevic, and D. Crabtree, MNRAS 313, 271 (2000). 

20. S.B. Pope, Turbulent flows (Cambridge University Press, 2000). 

21. N.N. Samus, E.V. Kazarovets, N.N. Kireeva, et al. General Catalogue of Variable Stars 
(2011). 

22. R. Stothers, Astrophys. J. 156, 541 (1969). 

23. R. Stothers, Astron. Astrophys. 18, 325 (1972). 

24. G. Wuchterl and M.U. Feuchtinger, Astron. Astrophys. 340, 419 (1998). 

25. M. Yang, B.W. Jiang, Astrophys. J. 727, 53 (2011). 

26. M. Yang, B.W. Jiang, Astrophys. J. 754, 35 (2012). 



10 



FIGURE CAPTIONS 



Fig. 1. Evolutionary tracks of red supergiants M Z ams — 2OM with metal abundances Z = 
0.004, 0.008 and 0.02 in the HR diagram. Parts of tracks corresponding to thermonuclear 
core helium burning are shown in solid lines, whereas parts of tracks with gravitational 
contraction as the energy source in the stellar center are shown in dotted lines. Arrows 
at the solid lines indicate the direction of evolution in the begining of helium burning. 
The mass loss factor / is given at the upper ending points of the tracks. 

Fig. 2. The kinetic energy growth rate rj as a function of the central helium abundance X( 4 He) c 
in hydrodynamical models Mzams = 15, 18, 20 and 25M with / = 0.3 and Z = 0.008. 
Hydrodynamical models are shown in filled circles. Models of each evolutionary sequence 
are connected by dotted lines. 

Fig. 3. Variations of the gas flow velocity at the upper boundary U (a) and bolometric magnitude 
M^oi (b) in the red supergiant with initial mass Mzams = 2OM , central helium abundance 
X( 4 He) c = 0.39 and pulsation period II = 383 days. 

Fig. 4. Same as Fig. 3 but for M ZAM s = 9M , X( 4 He) c = 0.36 and IT = 52 days. 

Fig. 5. Evolutionary changes of the luminoisity L and the period of radial pulsations II in the 
red supergiants Mzams = 2OM with metal abundances Z = 0.004, 0.008 and 0.02. 
Hydrodynamical models are represented by filled circles connected for each evolutionary 
sequence by solid lines for / = 0.3 and dotted lines for f — 1. Arrows indicate direction 
of evolutionary changes in the early stage of thermonuclear helium burning. 

Fig. 6. The period-luminosity diagram for red supergiants. Hydrodynamical models are rep- 
resented by filled circles [Z = 0.008), open circles [Z = 0.004) and open triangles 
[Z = 0.02). Linear fits ([6]) are shown by dashed lines. 

Fig. 7. The period-mass diagram for red supergiants with initial mass M Z ams = 2OM . Hy- 
drodynamical models are represented by filled circles connected by the dotted line for 
each evolutionary sequence [Z = 0.004 Z = 0.008). The arrows indicate the direction of 
evolution in the diagram. 

Fig. 8. The period-mass diagram for red supergiants of Magellanic Clouds. Hydrodynamical 
models with metal abandances Z = 0.008 are Z = 0.004 are shown by filled and open 
circles, respectively. Boundaries of the red supergiant strip fitted by relations (jSJ) are 
shown by dashed lines. 
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Figure 1: Evolutionary tracks of red supergiants M Z ams — 20M Q with metal abundances Z = 
0.004, 0.008 and 0.02 in the HR diagram. Parts of tracks corresponding to thermonuclear core 
helium burning are shown in solid lines, whereas parts of tracks with gravitational contraction 
as the energy source in the stellar center are shown in dotted lines. Arrows at the solid lines 
indicate the direction of evolution in the begining of helium burning. The mass loss factor / is 
given at the upper ending points of the tracks. 
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Figure 2: The kinetic energy growth rate r] as a function of the central helium abundance 
X( 4 He) c in hydro dynamical models M ZAM s = 15, 18, 20 and 25M Q with / = 0.3 and Z = 0.008. 
Hydrodynamical models are shown in filled circles. Models of each evolutionary sequence are 
connected by dotted lines. 
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Figure 3: Variations of the gas flow velocity at the upper boundary U (a) and bolometric 
magnitude Mbd (b) in the red supergiant with initial mass Mzams = 2OM , central helium 
abundance X( 4 He) c = 0.39 and pulsation period II = 383 days. 
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Figure 5: Evolutionary changes of the luminoisity L and the period of radial pulsations IT in 
the red supergiants M ZA ms = 20M Q with metal abundances Z = 0.004, 0.008 and 0.02. Hy- 
drodynamical models are represented by filled circles connected for each evolutionary sequence 
by solid lines for / = 0.3 and dotted lines for / = 1. Arrows indicate direction of evolutionary 
changes in the early stage of thermonuclear helium burning. 
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Figure 6: The period-luminosity diagram for red supergiants. Hydro dynamical models are 
represented by filled circles (Z = 0.008), open circles (Z = 0.004) and open triangles (Z = 0.02). 
Linear fits ()6]) are shown by dashed lines. 
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Figure 7: The period-mass diagram for red supergiants with initial mass Mzams — 2OM . 
Hydrodynamical models are represented by filled circles connected by the dotted line for each 
evolutionary sequence (Z = 0.004 Z = 0.008). The arrows indicate the direction of evolution 
in the diagram. 
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Figure 8: The period-mass diagram for red supergiants of Magellanic Clouds. Hydrodynamical 
models with metal abandances Z = 0.008 are Z = 0.004 are shown by filled and open circles, 
respectively. Boundaries of the red supergiant strip fitted by relations (IE]) are shown by dashed 
lines. 
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